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Принципы селективно-яркостной 
сегментации изменений 

внешнего облика объекта 
моноспектральной космической съемки 


В статье анализируются задачи, которые решаются на этапе обработки данных дистанционного зондирования 
Земли и предусматривают интерактивное участие человека в оценке качества материалов спутниковой 
съемки, улучшении полученных изображений и тематическом использовании в зависимости от 
разрабатываемых предложений. 


Введение 


Большое число задач, решаемых на этапе обработки данных дистанционного 
зондирования Земли, предполагает интерактивное участие человека в оценке качества 
материалов спутниковой съемки, улучшении полученных изображений и тематичес- 
ком использовании в зависимости от разрабатываемых приложений. 

Одной из таких задач является мониторинг поверхности Земли с обнаружением 
происходящих во времени изменений, вызванных теми или иными причинами. Тема- 
тическая обработка информации в таком применении заключается в сопоставлении 
разновременных изображений интересующего участка, обнаружении и классификации 
обнаруженных изменений. Существуют как традиционные (визуализация, фильтрация, 
получение новых изображений расчетными методами), так и современные методы 
такой обработки, использующие адаптивные алгоритмы на основе топографических 
отображений [1-3]. Последние применяются для более содержательной (с точки зрения 
тематической интерпретации) классификации и сегментации космических изображе- 
ний. Использование алгоритмов на основе топографических отображений позволяет 
добиться наилучших результатов при выделении тематической информации. 

В основе селективно-яркостной сегментации изменений, обнаруживаемых при 
сопоставлении разновременных изображений спутниковой съемки одного и того же 
наземного объекта, лежит анализ разностной гистограммы этих изображений [4]. Как 
известно, такие гистограммы представляются как следующая зависимость числа пик- 
селей от яркости: 


АН(Р) =Н(Р) - Но (Г), (1) 
где Н, (Е) — гистограмма изображения, полученная как характеристика «штатного» сос- 


тояния объекта, отклонения от которого в последующем и являются интересующими, 
представляемыми как изменения внешнего облика объекта; Н(Г) — гистограмма изоб- 


ражения текущей съемки объекта, рассчитанной на контроль возникновения возмож- 
ных отклонений его состояния от «штатного». 
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С учетом выполняемой роли гистограмма Ну, (Г) имеет название эталонной, опор- 


ной или базовой. 

Для реальных изображений такая гистограмма, соответствующая ряду измене- 
ний внешнего облика объекта за время между съемками, имеет вид, представленный 
на рис. 1. 
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Рисунок 1 — Разностная гистограмма разновременных изображений 
объекта спутниковой съемки 


Рис. 1 заставляет обратить внимание на то, что разностная гистограмма, при на- 
личии произошедших изменений внешнего облика объекта, состоит из группировок 
положительных и отрицательных значений площадей (группировок пикселей’) на 
отличающихся интервалах яркостей. Особенностями образования этих группировок 
являются изменения гистограммы текущего изображения по сравнению с гистограммой 
предшествующего ему эталонного изображения. В свою очередь упоминаемые измене- 
ния всегда оказываются следствием преобразований отражающих свойств компонент” 
рассматриваемого объекта на площадях тех компонент его эталонного состояния, ко- 
торые подверглись этим преобразованиям. 

Как уже было показано при оценке информационного содержания гистограмм [5], 
изменения отражающих свойств компоненты объекта естественно ведут к замене части 
или всей площади, занимаемой ею в составе этого объекта. Кроме того, такие изменения 
могуг охватывать частично или полностью площади не одной, а целого ряда компонент. 

В результате площадь объекта остается всегда неизменной. Однако материалы 
или элементы поверхности, которые формировали штатное распределение его яркостей 
на изображении, при изменениях внешнего облика объекта частично или полностью 
исчезают, уступая место другим. Поэтому наряду с положительными по знаку изме- 
нениями гистограмм, свойственными появлению новых по отражающим свойствам 
площадей компонент, появляются отрицательные изменения, отражающие уменьшение 
площадей компонент как следствие «оккупации» их территории новыми компонентами. 

При этом новым компонентам, как уже отмечалось, в природе и технике харак- 
терны распределения вероятностей отражающих свойств их поверхностей. А это ведет 
к тому, что изменения объектов носят не точечный, а интервальный характер, что и 
отражает разностная гистограмма. А границы интервала - это и есть те два порога, 
поиск которых открывал бы возможность яркостной сегментации изменений внешне- 
го облика объекта по его разновременным изображениям. 


' Следует иметь в виду, что каждому пикселю изображения, а также связанной с этим изображе- 
нием гистограммы, соответствует элемент земной поверхности, определяемый как пространственное 
разрешение (разрешающая способность) съемки. 

? Отражающие свойства компонент, составляющих объект, пропорциональны их яркостям на 
интересующих изображениях. 
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Принципы сегментирования 


Обратимся, в частности, к рассмотрению рис. 1. Здесь отметим, что положитель- 
ные группирования как раз и соответствуют здесь появлению в составе объекта за 
время между съемками новых компонент, характеризуемых яркостями каждого тако- 


го (р-го) интервала [^ „= Г) где р=1,2, 3, 4. Автоматическое (программное) 


ршш ртах › 
определение указанных интервалов представляет собой отдельную задачу, основанную 
на объединении в их составе яркостей группирования пикселей, соответствующих вы- 
полнению следующего условия: 


АН(Е)>0 УГ. 
Одновременному сокращению площадей старых компонент в общем случае со- 
ответствуют 4 интервалов разностной гистограммы [Ги = [вах › Содержащих отри- 


цательные группирования пикселей. Для той же ситуации изменений, представленной 
разностной гистограммой (рис. 1), имеем 4 =1. В результате, согласно изложенному 
следует считать, что за счет изменений одной компоненты объекта, имеющей яркос- 


ти также в легко определяемом диапазоне Гу, = [1 


Тит мах > 


здесь произошли яркостные 


изменения в четырех интервалах ( р =1, 2, 3, 4.) 

Таким образом, выделение яркостных интервалов изменений, возможность ко- 
торого открывает разностная гистограмма разновременных изображений наземного 
объекта, представляет собой первый принцип селективно-яркостного сегментирования 
изменений его внешнего облика. 

В целях дальнейшего анализа принципов использования разностных гистограмм 
обратимся к их моделированию. Причина обращения к моделированию состоит в от- 
сутствии точного априорного знания как пространственных характеристик изменений, 
так и распределений вероятностей яркостей изменяемых и измененных компонент 
для реальных изображений спутниковой съемки. Модели, предполагая известными 
интересующие характеристики, должны позволить зрительно и теоретически просле- 
дить особенности обнаружения изменений, представить их сущностные свойства, ко- 
торые внешне проявляются в разностных гистограммах экспериментальных съемок 
всего лишь как упомянутые выше группирования пикселей. 

С целью приближения такого моделирования к реальным условиям, позволяющим 
избежать предвзятости в формировании модели гистограммы, воспользуемся полу- 
ченным реально по данным спутниковой съемки односпектральным изображением 
участка земной поверхности, рассматриваемым как интересующий объект (рис. 2а). 
На рис 26 и в приведены соответственно гистограмма этого изображения, принятая в 
качестве эталонной Н.(Ё), и ее приближенная 6-ю компонентами модель Н(Г), иден- 
тифицированная с ней по приблизительному совпадению, что достаточно для обеспечения 
достоверности и наглядности результатов. Параметры моделирующих ограниченных 
нормальных распределений сведены в табл. рис. 2г. 

На полученной основе для последующего анализа предложена модель гистограммы 
текущей съемки того же объекта Н([) (рис. 29), представляющая изменения эталон- 
ной гистограммы путем включения при моделировании в ее состав 7-й компоненты — 


' Следует отметить, что в реальных условиях определение интервала группирований до пиксе- 
ля АН(Г) часто затруднено в силу влияния шумов изображений, вызывающих существенное для этой 
задачи пропадание пикселей. Это явление особенно чувствительно для небольших значений на грани- 
цах интересующего интервала АН(Ё„) или АН(Ё„). 


шш 
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компоненты изменений. Для удобства анализа принято, что она имитирует результат 
1%-го замещения площади одной из компонент (№ 6) эталонного изображения, зани- 
мающей до этого 14% площади объекта. В качестве параметров ее распределения 
выбраны х, =55 и о =4. Разностная гистограмма моделей текущей и эталонной гисто- 


грамм показана на рис. 2е. 
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Рисунок 2 — Исходные данные и результат моделирования разностной гистограммы 
(а — изображение участка космической съемки размером 128898 пикселей, 
б — его гистограмма, в — модель эталонной гистограммы, г — таблица компонент параметров 
распределений модели эталонной гистограммы, д -— модель текущей гистограммы 
с изменениями, е — разностная гистограмма) 


Одновременно с этими результатами на рис. 3 даны взвешенные для суперпо- 
зиции распределения яркостей всех компонент, параметры которых и обеспечили при- 
близительную идентификацию модели, приведенной на рис. 28. 

Таким образом, сама разностная гистограмма обеспечивает, как это следует из 


> ЗЕ . Е 
рассмотрения рис. 2е, определение в составе объекта интервала яркостей Гути = ГСутах 


новой компоненты изменений (р =1). Одновременно рис. 3 позволяет установить, что 


распределение этой компоненты в пределах того же интервала яркостей пересекает- 
ся, по крайней мере зрительно, с распределениями 2-й, 3-й и 4-й компонент. 
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При этом здесь легко представить, что степень такого пересечения тем больше, 
чем больше дисперсии пересекающихся распределений и чем больше относительные 
площади, которые занимают описываемые компоненты. 


| [компонента 1 Ш —— компонента 1 


—— компонента 2 


число пикселей 


— компонента изменений 


—— компонента 5 
компонента 8 


1 29 57 85 13 141 163 197 225 253 
градациияриести 
Рисунок 3 — Набор распределений вероятностей яркостей компонент, 
взвешенных их относительной площадью в составе объекта и приведенных к числу пикселей 


его изображения (1-6 — распределения компонент штатного состава объекта; 
между 2-м и 3-м — распределение компоненты изменений) 


В этих условиях можно теперь обратить внимание на то, что выделенному по раз- 
ностной гистограмме интервалу яркостей соответствует площадь компоненты изме- 
нений: 

[= них 
5м = 5: У`АН(Е), 6) 
[= п пи 
где Ду — пространственная разрешающая способность спутниковой съемки; 
1. +. =43+69 — интервал изменений, определенный по разностной гисто- 


там тах 
грамме (положительных значений); ДН(Г.) — разностная гистограмма. 

С другой стороны, в силу особенностей моделирования гистограмм, истинным 
значением площади компоненты изменений должна быть следующая величина, опре- 
деляемая непосредственно по распределению вероятностей самой компоненты: 


изм = 4. Уз | ие) < ря ый (3) 
УЁ АН(Ё)>0 


где /,„(Ё)- плотность распределения т. яркостей компоненты изменений; 


5 изм 


— относительная площадь компоненты изменений (5 — площадь компоненты, 


изм 


5 — площадь объекта); М — число пикселей изображения объекта; [...] — означает це- 
лую часть рассматриваемого выражения. 
Отсюда, принимая во внимание сноску на странице 2, оценку площади изменений 


на интервале [лиш = [ль ‚ следует считать приближенной: 
Ех 5 
Зизм =Ду. уз Ё к Чан. |- Ду. уз о-ва (4) 
Е, АНЕ>0 РИ 5 
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В то же время при моделировании без имитации шумов изображения это равенство 
оказывается точным, так как при определении границы интервала ни один пиксель 
группирования изменений не исчезает. В частности, преобразовав выражение (2) в 
соответствии с принципами моделирования гистограмм [6], получим: 


ВЕБ Е = пах ВЕЕТ ная 
би = Д5- УАН(Е) = 4х. УН(Е)- УНЬ (Е) = 
= па Е = п = па 


= 5. > У лю |- у РИО = 2 


= вы Е #1 = ь В 7 


=бае ОН 


1=Г пы 


Как раз в последнем выражении (5) при сопоставлении уменьшаемого и вычи- 
таемого аналитически наглядно прослеживается компенсация «хвостов» распределений, 
соседствующих с распределением компоненты изменений и выделенных изучаемым 
интервалом яркостей. 

Полученный результат — легко достигаемый уровень сегментации изменений 
внешнего облика объекта. В то же время знания интервала яркостей и размеров (площа- 
дей) сегментов изменений без определения их мест на изображении, а следовательно, 
на территории объекта обычно недостаточны для решения задачи обнаружения из- 
менений. 

Для выбора пути устранения указанного недостатка необходимо обратиться к 
уточнению особенностей связи в паре гистограмма-изображение. Она достаточно оче- 
видна. Если гистограмма является яркостным распределением пикселей изображения 
НС), то само изображение представляет собой их пространственно-яркостное распре- 


деление [(х, у). Поэтому именно яркостная связь этих двух распределений открывает 


возможность определения пространственных положений (координат) элементов разре- 
шения каждого найденного сегмента изменений. 


= 
В частности, применив выделение яркостей пикселей р-го интервала Г, С =, м 


на текущем изображении, удается легко устранить упомянутый информационный 
недостаток и получить знания о том, как распределены в пространстве изображения 
эти яркости после изменений объекта. Причем наиболее целесообразно это преобра- 
зование выполнять путем бинаризации и представления результата в виде множества: 


1 Их, ут -р . 
В (х = при (х, У ) «(2 ‘ри ти (6) 
.. 0 при Цх, У,0) = (р Я ты ы 


Возвращаясь вновь к рассматриваемому примеру моделирования разностной 
гистограммы, открываем, что выделение частей текущего изображения, имеющих яр- 


кости в единственном здесь интервале изменений у я Уре ‚ должно иметь итогом 


Тит 


представление их в виде 7-ми составляющих сегментов, так как согласно (5) соот- 
ветствующая площадь в конечном итоге определяется следующим образом: 


каю УВЕА ЫР ‚м -& [и РО Ум 


= ур 1 о 1 Ё ПА 
РЕН них РЕБ ах р. (7) 
+45. У. [ла м. +45. У) |1 ы 
[= о [= т 
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Очевидно, что каждая из получаемых здесь составляющих площадей представ- 
ляет собой часть компонент объекта, имеющих яркости в рассматриваемом диапазоне. 
Все они образуют множество, состоящее из разнесенных (в общем случае) в пространстве 
измерений отдельных подмножеств. Каждое из них должно иметь свое геометрическое 
положение в составе соответствующего изображения. Отсюда текущее изображение 
объекта, преобразованное согласно (6), для особенностей приведенных моделей гисто- 
грамм, должно иметь вид примерно такой, как показан на рис. 4. 


Рисунок 4 — Типовая картина бинаризованного текущего изображения 
в заданном яркостном интервале 


Известно, что гистограммы спутниковых изображений наземных объектов в 
большинстве случаев представляют собой распределения площадей их компонент по 
яркостям, не сосредоточенные на одном избранном интервале. Поэтому яркостные 
характеристики компонент изменений чаще всего могут совпадать с яркостными ха- 
рактеристиками некоторых компонент объекта. Такая, в частности, ситуация характерна 
для моделируемых гистограмм (рис. 2в, д), где согласно рис. 3 части распределений ве- 


роятностей компонент («хвосты») имеют яркости, которые входят в интервал [лиш = [1 


Тит Т тах * 

Именно представляемые этими «хвостами» части компонент (каждой компоненте в 

целом соответствует все распределение) наряду с сегментом изменений составляют 
бинаризованное изображение рис. 4. 

Конечно, если не осуществлять введенной при этом бинаризации, то каждый из 

сегментов текущего изображения вместо однотонных и максимально ярких здесь их 

представлений мог бы иметь отличающиеся яркости на черном фоне. Но все они обра- 


зовывали бы вместе выделяющий их и заранее известный интервал лиш = 1х = 43+ 69. 


Новой информации это не позволило бы получить. Поэтому, интересуясь положением 
этих сегментов в пространстве изображения, от яркостей можно отказаться. Одновре- 
менно, переход к бинаризации в этой ситуации позволяет, во-первых, избежать плохой 
наблюдаемости рассматриваемых сегментов при невысоких яркостях интервала. Во- 
вторых, выделенные при этом сегменты — это составляющие единого множества, 
операции, которые достаточно удобны в последующих преобразованиях. 

Таким образом, вторым принципом локально-яркостной сегментации является 
переход от яркостных представлений об изменениях (от разностной гистограммы) 
к пространственным представлениям о размещении всех сегментов текущего изобра- 
жения, характеризующих в интервале сегментации состояние объекта после изме- 
нений его внешнего облика. 

Здесь следует обратить внимание на то, что рассмотренное до этого определение 
яркостного интервала изменений и размеров участка, на котором они произошли, не 
зависит от геометрических несоответствий исходных изображений друг другу, так 
как гистограммы не связаны с пространственными характеристиками изображений. 
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В отличие от этого любые искажения геометрической передачи пропорций и размеще- 
ния компонент объекта, свойственные текущему изображению, отражаются на раз- 
мещении компонент и их частей, характеризуемых яркостями изучаемого интервала 
([^. -Г..). Однако при отсутствии сведений о геометрическом смещении или повороте 


рт ртах 

текущего изображения упомянутые искажения внешне никак здесь себя не проявляют. 
К сожалению, выполненное пространственно-яркостное выделение на текущем 

изображении еще не позволяет отделить сегменты изменений от всех наблюдаемых 

здесь сегментов заданного интервала (рис. 4). Для устранения такого недостатка 

достаточно обратиться к эталонному изображению. Применение к нему выделения ком- 

понент и их частей, аналогичного реализованному для текущего изображения, пред- 

ставляется следующей бинаризацией: 


В°(х,у)= при 7й о(х, уд ЕЕ р ит ПЕ (8) 
0 при (ху, (Е а 

В соответствии с принятой к рассмотрению моделью гистограммы эталонного 
изображения число компонент объекта в заданном яркостном интервале составляет 6 
в отличие от текущего изображения. Эталон характеризует штатное состояние объекта 
до любых не запланированных его изменений. Поэтому здесь компонента изменений 
отсутствует. В итоге суммарная площадь всех имеющих место компонент оказывает- 
ся равной: 


45. У Н.(Р)=4у. У. 


1-х 1-х р 


1=Н тах 
Е ум = Ду. > [ме м 


И И, = Я Е Е (9) 
ТЕГ их ТЕ тах 
+45. У. ла 2 № +. +5. У. [дах ы 
= ты 1-Е ты 


Согласно этому преобразованное эталонное изображение должно выглядеть сле- 
дующим образом (рис. 5). 


Рисунок 5 — Типовая картина бинаризованного эталонного изображения 
в заданном яркостном интервале 


Отсюда третий принцип локально-яркостной сегментации состоит в переходе 
от яркостных представлений об изменениях (разностной гистограммы) к пространст- 
венным представлениям о размещении всех компонент эталонного изображения, 
имеющих яркости в заданном интервале и отражающих состояние объекта до из- 
менений его внешнего облика. 

Точно так же, как и для текущего изображения, здесь любые искажения геомет- 
рической передачи пропорций и размещения компонент объекта, свойственные эта- 


лонному изображению, отражаются в размещении обнаруживаемых компонент и их 
частей, характеризуемых яркостями интервала изменений ( [7 „+1 ). Отличия здесь 


шт тах 
р 


состоят только в том, что это изображение получено в другое время. Поэтому геомет- 
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рические искажения его должны, как правило, не совпадать с искажениями текущего 
изображения. В то же время при рассмотрении изолированно только этого результата 
бинаризации упомянутые искажения себя и здесь не проявляют. 

Логически следующим этапом локально-яркостной сегментации является сопо- 
ставление яркостно бинаризованного текущего изображения с эталонным с целью 
определения, какие сегменты этого интервала сохранились на нем за время между 
съемками и какие оказались приобретенными. Аналитически выполняемой для этого 
операцией является разность множеств (6) и (8), состоящих из единичных по яркостям 
пикселей, каждый из которых характеризуется координатами положения в пространстве 
изображения: 


В, (х,у) = В,(х,у)\ Вь(х,у). (10) 


Рисунок 6 — Типовая картина разности бинаризованных текущего 
и эталонного изображений в яркостном интервале обнаружения изменений 


Наряду с этим, выполнению этой операции соответствует вычитание изображе- 
ний согласно следующему правилу: 


1 при В, (х, у) =1л Вь(х, у) = В 


Вл (х,у) = В,(х,у)-Вь(х,у) = (11) 


0 при всех остальных сочетаниях значений. 


Указанная разность при отсутствии взаимных геометрических искажений текущего 
и эталонного изображений сегментов рассматриваемого яркостного интервала (рис. 4 и 5) 
приводит к выделению сегмента изменений и имеет следующее представление (рис. 6). 

Таким образом, четвертый принцип локально-яркостной сегментации состоит 
в определении пространственного положения сегмента изменений на основе исполь- 
зования разности бинаризованных изображений в яркостном интервале обнаруженных 
изменений (1 -+Ё,)- 


риш ® Грпах 
Ограничениями таких возможностей здесь и являются часто имеющие место гео- 
метрические искажения эталонного и текущего изображений. Если же они одинако- 
вы для того и другого, то рассмотренный результат будет характеризоваться только 
ошибками пространственного положения компоненты изменений, демонстрируемой 
рис. 6. Здесь они отражались бы в некотором смещении в плоскости изображения одно- 
го единственного сегмента — области изменений. Однако, если наличие или отсутствие 
искажений нечем подтвердить, то зрительно это не может быть никак обнаружено. 
При отличиях геометрических искажений пары сопоставляемых изображений, 
кроме смещенности положения сегмента изменений, в разности их возникают новые 
сегменты, образуемые остатками от рассмотренного выше вычитания компонент би- 
наризованных изображений, не подвергнутых изменениям, но не совпадающих теперь 
друг с другом. Так, на рис. 7 приведены эталонное и текущее изображения, бинаризо- 
ванные в интервале яркостей при имеющемся взаимном геометрическом несоответствии 
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(текущее изображение смещено влево по сравнению с эталонным) и приблизительный 
вид их разности. 


Рисунок 7 — Геометрически смещенное текущее изображение объекта (а), 
эталонное изображение (6) и их разность (в) 


К аналогичным разностям приводят любые геометрические искажения сопостав- 
ляемых изображений, бинаризованных в яркостном интервале обнаруженных измене- 
ний. И чем больше геометрические сдвиги разновременных изображений, тем большие 
площади приобретают получаемые здесь разностные остатки. Однако сегмент истинных 
изменений при этом сохраняет свое присутствие в составе разностного изображения. 
Теперь он, характеризуясь смещенными координатами в силу тех же искажений, тре- 
бует еще и своего обнаружения среди указанных остатков. Последние на разностном 
бинарном изображении (рис. 7в), выступают как отдельные сегменты изменений в 
рассматриваемом яркостном интервале. Пространственное выделение каждого из 
них играет роль ложной тревоги в оценке изменений внешнего облика объекта. 

Как уже отмечалось в начале рассмотрения принципов локально-яркостной сег- 
ментации, в составе разностной гистограммы (рис. 2е) наряду с положительными 
группированиями пикселей присутствуют отрицательные группирования, отражающие 
уменьшение площадей компонент как следствие «оккупации» их территории новы- 
ми компонентами. 

Если вначале интервал яркостей, соответствующий отрицательным значениям 
разностной гистограммы, не представлял интереса, то теперь следует обратить на 
него внимание. В рассматриваемой ситуации моделирования этот интервал равен 
Вы и = 90255. 

Учитывая простоту ситуации, обусловленную всего одним изменением объекта 
в результате 1%-го замещения площади всего одной из его компонент, и руководст- 
вуясь особенностями разностной гистограммы, достаточно просто увидеть, что ука- 
занный интервал яркостей характеризует исходные яркости конкретной компоненты 
эталонного изображения. Именно эта компонента подверглась изменениям, обнаружи- 
ваемым на текущем изображении по разностной гистограмме. 


404 «Искусственный интеллект» 472009 


Принципы селективно-яркостной сегментации изменений внешнего облика... 8 Б 


Отсюда, если выполнить бинаризацию эталонного изображения в интервале 


исходных яркостей изменяемой компоненты ( ), то результатом ее должно 


а тт + 4 тах 
быть изображение пространственного положения этой компоненты. 

Бинаризация здесь выполняется по правилам, уже известным из предыдущего 
изложения: 


1 при У (х, ув (Е тт гро: (12) 
0 при [4 (х, д, =. 


шт 4 тах 


0 
В. (х,у)= 


Результат бинаризации в диапазоне им = 1х =90+255 в принципе должен иметь 
представление, показанное на рис. 8. 


Рисунок 8 — Бинаризованное изображение компоненты, 
подвергшейся изменениям 


Часть площади рассматриваемой компоненты (1%) согласно предложенным усло- 
виям моделирования должна занять компонента изменений с яркостью в диапазоне 
И + А 


Тпит пах * 


Ее бинаризованное представление в том же пространстве уже приведено 


на рис. 6 и в составе компонент рис. 7в. 
РИ 


| тах 


Таким образом, два отличающихся между собой по яркостям интервала //, 


= 


пах 


Пит 


и Ё 


т характеризуют с разных сторон одну и ту же пространственно располо- 
женную компоненту. До того, как изменения произошли, вся площадь, которую за- 
нимала эта компонента, имела одни яркости, а затем часть ее изменила их, но при 
этом сама компонента не изменила своего положения. Это дает основания для проверки 
совместимости изображений (компоненты изменений и компоненты, подвергнутой 
изменениям), полученных двумя разными путями с помощью операций пресечения 
соответствующих множеств: 


1 придВ, (х, у) =1^ В (х, у) =1; 
Ва(х у) = ВА (х, УП В (с, у)= - : 


при всех остальных сочетаниях значений. 


(13) 


При этом вполне очевидно, что значимый результат (непустое множество) от 
пересечения рассматриваемых множеств, представленных изображениями рис. 6 и 8, 
однозначно свидетельствует, что компонента изменений (рис. 6) появилась за время 
между съемками на площади компоненты эталона, изображенной на рис. 8, наглядным 


и ожидаемым результатом представления множества В; (х, у) в этом случае является 


изображение, полностью совпадающее с изображением рис. 6. 


'Рассмотренное представление не отражает соотношения площадей изменяемой компоненты 
объекта и компоненты изменений, а преследует целью только проиллюстрировать выполняемые со- 
поставления. 
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Такой итог обработки достаточно очевиден при отсутствии взаимных геомет- 
рических искажений пространства изображений. Не принимая этого во внимание, от 
операции пересечения множеств можно отказаться: она не внесла ничего нового в ре- 
зультат обработки, так как 

Ва (х, у) = ВА (х, у). (14) 

Однако, если иметь дело с анализируемыми выше геометрическими искажения- 
ми изображений (текущего и эталонного), то ситуация существенно меняется. Вместо 
изображения изменений (рис. 6) приходится располагать таким изображением, но с 
помехами (рис. 7в). 

Характерно, что пересечение множества, представленного таким изображением, 
содержащим бинаризованные сегменты остатков как результат геометрических иска- 
жений (рис. 76), с множеством изображения (рис. 8) приводит практически к тому же 
результату, то есть к изображению рис. 6. 

Причина идентичности здесь состоит в том, что изображение (рис. 8) не может 
содержать каких-либо компонент на местах сегментов остатков, свойственных изобра- 
жению рис. 7в. Если рис. 7в образован как разность В, (х, у) геометрически несведен- 


ных текущего В,(х, у) и эталонного В, (х, у) изображений, то рис. 8 представляет только 


[ 
4-е яркостное выделение В, (х, у) из эталонного изображения. 


К сожалению, следует заметить, что те же геометрические искажения могут и здесь 
влиять на результат В; (х, у). Так, множество обнаруженного сегмента изменений В) (х, у) 


0 
в силу его смещения может быть только частью множества В й (х, у). Несовпадающая 


часть его при пересечении будет потеряна, и может сложиться представление о мень- 
шей степени изменений, чем они имеют место фактически. В то же время для пропуска 
изменений необходимо иметь геометрическое смещение, приводящее к несовпадению 
рассматриваемых множеств (пустое множество от пересечения). 

Одновременно можно обнаружить, что при расположении замещаемой части ком- 
поненты изменений внутри компоненты, подвергаемой изменениям, какое-либо умень- 
шение размеров вследствие геометрических искажений может не проявиться вообще. 

Отсюда пятый принцип локально-яркостной сегментации состоит в избавле- 
нии от ложных тревог сегментации, вызываемых геометрическими искажениями 
разновременных изображений, путем образования пересечения множества измене- 
ний с множеством локально-яркостного представления и бинаризации эталонного 
изображения в интервале яркостей компоненты, подвергнутой изменениям. 

Комментируя этот принцип локально-яркостной сегментации можно отметить, что, 
кроме образования помех при геометрических различиях разновременных изображений, 
на изображениях, которые участвуют в формировании множеств, присутствуют шумовые 
яркостные выбросы. Они обусловлены физическими свойствами фотосенсоров матрицы 
спутникового устройства регистрации светового потока, отраженного поверхностью 
Земли, а также шумами электронных компонент и трактов передачи и наземного приема 
изображений. При этом на разновременных изображениях (эталонном и текущем) они 
независимы. Отсюда, возвращаясь к операции пятого принципа локально-яркостной сег- 
ментации (11), обнаружим, что разность В‚(х, у) оставляет здесь нескомпенсированны- 


ми шумы выбросы яркости только текущего изображения В, (х, у). Реализация очеред- 
ного, шестого, принципа В+ (х,у) = В (ху) ПВ» (х, У), то есть образование пересечения в 
силу принадлежности шумов изображения В? (х, У) эталонному (более раннему по вре- 


мени) изображению удаляет шумовые выбросы разности В„(х, У) в силу несовпадения 
мест их проявления на сопоставляемых при этом изображениях. 
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Выполняемое пересечение бинаризованного разностного изображения с бина- 
ризованным изображением эталона в диапазоне компонент, подвергаемых измене- 
ниям, обеспечивает не только избавление от мешающих факторов геометрической 
некорректности, но и от шумовых выбросов, свойственных космической съемке. 


Экспериментальная демонстрация применения 


Для демонстрации результатов практического применения принципов локаль- 
но-яркостной сегментации на рис. 9 представлены односпектральные изображения 
спутниковой съемки Гап@зай 7 участка поверхности Земли (Япония) до и после зем- 
летрясения и результаты их обработки на каждом из двух интервалов обнаруживае- 
мых изменений. 


а) эталонное изображение (съемка 13.08.04) 6) текущее изображение (съемка 24.10.04) 


число пикселеи 


градации яркости : : : 
г) бинаризованный эталон 


в) разностная гистограмма изображений в интервале 36-56 ( В°(х.5))) 
а > 


д) разность В», (х, у) в интервале 13-33 е) пересечение изображений В, (х, у) П В, (х, у) 


ж) разность В„› (х, у) в интервале140-152 3) пересечение изображений В,, (х, у) Г] В, (х, у) 


Рисунок 9 — Демонстрационная обработка реальных изображений спутниковой съемки 


Выводы. Первый локальный яркостный интервал изменений (Д„-=Д шк =13-33) 
(5 = 3254) оказался соответствующим в основном шумам или ошибкам геометри- 
ческого сведения, которые при образовании разности В), (х, у) еще проявлялись (рис. 90) 
и большей частью были скомпенсированы после пересечения изображений (рис. 9е). 


Второй локальный яркостный интервал изменений ( [2 = [2х =140-152) (з= 136253) 
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оказался информативным, позволяющим после компенсации шумов (рис. 93) обнару- 
жить, что выделенным сегментом соответствующих изменений является вся река; 
при этом разностная гистограмма (рис. 9в) свидетельствует о том, что яркость этого 


сегмента увеличилась по сравнению с исходным ее значением Ги = Дик = =36-— 


Тит 1 тах 
— 56 (5=156818, а должно быть 168794) до яркостей [1 +[„ =140-152, что мо- 


жет означать возрастание коэффициента отражения воды вероятнее всего в результате 
подъема со дна реки иловых отложений и размытия берегов при землетрясении. 


Заключение 


Проведенное рассмотрение позволяет сделать следующие выводы: 

1. Локально-яркостная сегментация изменений внешнего облика объекта косми- 
ческой съемки строится на утверждении: если разностная гистограмма разновременных 
изображений объекта космической съемки ДН (Г) = Н([)-Н.(Ё) имеет положительные 


значения в р-ом интервале яркостей [/ „= то в том же порядке образованная 


ршш ришах ?› 


разность самих изображений [(х,у) и Г, (х, у), предварительно бинаризованных в 
> 0 
этом же (р-ом) интервале яркостей В, (х, у), В,(х, у) ‚ и рассматриваемая как множест- 


во В\(х, у) =В,(х, у) \ В (х, у), представляет собой отображаемую область объекта 


(сегмент) в виде множества значений координат пикселей этих яркостных изменений 
за время между сеансами съемки. 

2. Пять принципов сегментации, пошагово реализуемые в процессе достижения 
конечного результата, позволяют определить геометрическое расположение искомых 
изменений, диапазон их яркостей, избавиться от помех, вызывающих ошибки лож- 
ной тревоги, и защититься от недостаточной точности геометрического сведения изо- 
бражений разновременной космической съемки. 

3. Выполненное экспериментальное рассмотрение применения принципов локаль- 
но-яркостной сегментации изменений дает возможность в приведенной практичес- 
кой ситуации увидеть и оценить существо каждого шага их осуществления. 
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Л.А. Блозерський, Л.В. Орешюна 

Принципи селективно-яскрав!снот сегментацй зм!н зовнпинього виду 0б’екта 

моноспектрально! косм1чно! зйомки 

У стати! розглядаються задач, як! виршгуються на етатш обробки даних дистаншйного зондування Земль, 
та передбачають 1нтерактивну участь людини в ощнц! якост! матер1алв супутниково! зйомки, покращенн! 
отриманих зображень та тематичному використанн! в залежност! в1д думок, що розроблюються. 
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